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Abstract. We continue the study of the Hermitian random matrix ensemble with external 
source 

where A has two distinct eigenvalues ±a of equal multiplicity. This model exhibits a phase 
transition for the value a = 1, since the eigenvalues of M accumulate on two intervals for 
a > 1, and on one interval for < a < 1. The case a > 1 was treated in part I, where it 
was proved that local eigenvalue correlations have the universal limiting behavior which is 
known for unitarily invariant random matrices, that is, limiting eigenvalue correlations are 
expressed in terms of the sine kernel in the bulk of the spectrum, and in terms of the Airy 
kernel at the edge. In this paper we establish the same results for the case < a < 1. As in 
part I we apply the Deift/Zhou steepest descent analysis to a 3 x 3-matrix Riemann-Hilbert 
problem. Due to the different structure of an underlying Riemann surface, the analysis 
includes an additional step involving a global opening of lenses, which is a new phenomenon 
in the steepest descent analysis of Riemann-Hilbert problems. 



1. Introduction 

This paper is a continuation of [Hj to which we will frequently refer in this paper. It will 
be followed by a third part [Hj, which deals with the critical case. In these papers, we study 
the random matrix ensemble with external source A 

^^(^dM) = — e-^^'^^^'^^-^^'UM, (1.1) 

Zn 

defined on n x n Hermitian matrices M, with Gaussian potential 

V{M) = ^M^ (1.2) 

and with external source 

A = diag(a, . . . , a, —a, . . . , —a). (1.3) 

n/2 n/2 

In the physics literature, the ensemble (jl.lll was studied in a series of papers of Brezin 
and Hikami P-jin], and P. Zinn- Justin [S21, [SSI- Our aim is to obtain rigorous results on 
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eigenvalue correlations in the large n-limit using the steepest descent / stationary phase 
method for Riemann-Hilbert (RH) problems jJH]; thereby extending the works 0; ll], \^^ ■ 
[T7j . fIE\ . j2E] who treated the unitary invariant case (i.e., A = 0) with RH techniques. 

While the unitary invariant case is connected with orthogonal polynomials |T51I2Z|; the 
ensemble is connected with multiple orthogonal polynomials 0. These are character- 
ized by a matrix RH problem , and the eigenvalue correlation kernel of has a direct 
expression in terms of the solution of this RH problem, see [H], ^1] and also formula ()1.14j) 
below. The RH problem for (jl.lj) has size (r + 1) x (r + 1) if r is the number of distinct 
eigenvalues of A. So with the choice (jl.3p . the RH problem is 3 x 3- matrix valued. 

The asymptotic analysis of RH problems has been mostly restricted to the 2x2 case. 
The analysis of larger size RH problems presents some novel technical features as already 
demonstrated in [B] and [23 • the present paper another new feature appears, namely at 
a critical stage in the analysis we perform a global opening of lenses. This global opening of 
lenses requires a global understanding of an associated Riemann surface, which is explicitly 
known for the Gaussian case p.2|) . This is why we restrict ourselves to ()1.2j) although in 
principle our methods are applicable to more general polynomial V. 

The Gaussian case has some special relevance in its own right as well. Indeed, first of all 
we note that for (|1.2|) we can complete the square in and then it follows that 

M = Mo + A (1.4) 

where Mq is a GUE matrix. So in the Gaussian case the ensemble ()1.1|) is an example of a 
random + deterministic model, see also ^-jl3j. 

A second interpretation of the Gaussian model comes from non-intersecting Brownian 
paths. This can be seen from the joint probability density for the eigenvalues of M, which 
by the HarishChandra/ltzykson-Zuber formula ^HI, takes the form 

1 " 

J- n (^.-^fc)det(e"^^-);^^^^n^"^"'^' (1-5) 

" l<j<k<n j=l 

for the case ()1.2|) . Here ai, . . . ,a„ are the eigenvalues of A, which are assumed to be all 
distinct in ()1.5|) . In the case of coinciding eigenvalues of A we have to take the appropriate 
limit of ()1.5|) . see formula (3.17) in [S]. 

Formula p.5|) also arises as the distribution of non-intersecting Brownian paths. Consider 
n independent Brownian motions (in fact Brownian bridges) on the line, starting at some 
fixed points Si < S2 < ■ ■ ■ < s„ at time t = 0, ending at some fixed points 6i < 62 < . . . < 
at time t = 1, and conditioned not to intersect for t G (0, 1). Then by a theorem of Karlin 
and McGregor j2Tj, the joint probability density of the positions of the Brownian bridges at 
time t G (0, 1) is given by 

Pn{xi, ...,Xn) = -^ det{p{sj, Xk, t))lf,^^ det{p{xj, bk; 1 - t))lk=i (1-6) 

where p{x,y;t) is the transition kernel of the Brownian motion and C„ is a normalization 
constant. Let us consider a scaled Brownian motion for which 
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and let us take a limit when all initial points Sj converge to the origin. In this case formula 
fll.6|) takes the form 



Pnixi,...,Xn) = [[ {xj - Xk) det ie— \ ||e"5J(T^^^. (1.^ 
This coincides with (|1.5|) if we make the identifications 



A,- = — , , ak = bk\l — - — . (1.9) 

So at any time t G (0, 1) the positions of n non-intersecting Brownian bridges starting at 
and ending at specified points are distributed as the eigenvalues of a Gaussian random 
matrix with external source. 

The connection between random matrices and non-intersecting random paths is actually 
well-known, see e.g. the recent works j2], [201; |2H] and references cited therein. 

P. Zinn- Justin showed that the m-point correlation functions for the eigenvalues of M 
have determinantal form 

RmiXu . . . , A„) = det (ir„(A„ Xk)),<j^k<n. ■ (1-10) 

It was shown in j3] that the average characteristic polynomial 

P{z) =E[det{zI - M)] 

is a multiple orthogonal polynomial of type II, which for the Gaussian case ()1.2|) is a multiple 
Hermite polynomial, see PP, [7], and that the correlation kernel Kn can be expressed in terms 
of the solution of the RH problem for multiple orthogonal polynomials ^T], see also [T^ . 

We state the RH problem here for the Gaussian case (jl.2p and for the external source (jl.3j) 
where n is even. Then the RH problem asks for a 3 x 3 matrix valued function Y satisfying 
the following. 

• F : C \ R ^ C^^Ms analytic. 

• For X G M, there is a jump 

(1 wi{x) W2{xy 

1 I (1.11) 

1 



where 



and Y^{x) {Y-{x)) denotes the limit of Y{z) as 2; ^ a; G M from the upper (lower) 

half-plane. 

As z — >• cxD, we have 

\ 

Y{z) = {I + 0{l/z)) \ . (1.13) 

^-"/V 
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The RH problem has a unique solution in terms of multiple Hermite polynomials and their 
Cauchy transforms [HI]. The correlation kernel Kn is expressed in terms of Y as follows 

Kn{x,y) = ^ (0 e-'^y e""'^^) Y-\y)Y{x) I I (1.14) 

Our goal is to analyze the above RH problem in the large n limit and to obtain from 
this scaling limits of the kernel ()1.14j) in various regimes. In this paper we consider the case 
< a < 1. The case a > 1 was considered in and the critical case a = 1 will be considered 
in [Hj. First we describe the limiting mean density of eigenvalues. 

Theorem 1.1. The limiting mean density of eigenvalues 

p{x) = lim —Kn{x,x) (1-15) 

n^oo n 

exists for every a > 0. It satisfies 

p(x) = - |Im^(x)| , (1.16) 

71 

where ^ = C,{x) is a solution of the cubic equation, 

- x^^ - {a^ - 1)^ + xa^ = 0. (1.17) 

The support of p consists of those x G M for which (|1.17|) has a non-real solution. 

(a) ForO < a < 1, the support of p consists of one interval [—zi, zi], and p is real analytic 
and positive on {—zi,zi), and it vanishes like a square root at the edge points ±zi, 
i.e., there exists a constant pi > such that 

p(x) = — Ix =F 2i|^/^(l + o(l)) as X ^ ±zi, X e {-zi,zi). (1.18) 

TT 

(b) For a = 1, the support of p consists of one interval [—Zi,Zi], and p is real analytic 
and positive on {—Zi, 0) U (0, Zi), it vanishes like a square root at the edge points ±zi, 
and it vanishes like a third root at 0, i.e., there exists a constant c > such that 

p(a;) =c|x|^/^(l + o(l)), asx-^0. (1.19) 

(c) For a > 1, the support of p consists of two disjoint intervals \—z\, —Z2\ U [z2,zi\ with 
< Z2 < Zi, p is real analytic and positive on (— Zi, —-22) U {z2,zi), and it vanishes 
like a square root at the edge points ±zi, ±Z2. 



Remark: Theorem 11.11 is a very special case of a theorem of Pastur on the eigenvalues 
of a matrix M = Mq + A where Mq is random and A is deterministic as in ()1.4j) . Since 
in this paper our interest is in the case < a < 1, we show in Section 9 how Theorem 11.11 
follows from our methods for this case. See 161 for the case a > 1. 



Remark: Theorem 1 1.1 1 has the following interpretation in terms of non- intersecting Brownian 
motions starting at and ending at some specified points bj. We suppose n is even, and 
we let half of the fo/s coincide with 6 > and the other half with —b. Then as explained 
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before, at time t G (0, 1) the (rescaled) positions of the Brownian paths coincide with the 
eigenvalues of the Gaussian random matrix with external source ()1.3|) where 



^ (1.20) 



The phase transition at a = 1 corresponds to 

So, by Theorem ll.il the limiting distribution of the Brownian paths as n oo is supported 
by one interval when t < and by two intervals when t > t^. At the critical time tc the two 
groups of Brownian paths split, with one group ending at t = 1 at 6 and the other at —b. 

As in |6;^ we formulate our main result in terms of a rescaled version of the kernel Kn 

kn{x, y) = e'^('^(^)-'^(^))if„(x, y) (1.21) 

for some function h. The rescaling 1)1.211) does not affect the correlation functions p.lO|) . 

Theorem 1.2. Let < a < 1 and let z\ and p be as in Theorem (a). Then there is a 
function h such that the following hold for the rescaled kernel f)1.21|) . 

(a) For every Xq G (— -Zi, Zi) and u,v eM., we have 

1 - / u V \ sm7i(u — v) 

hm ^--Kn{xo + ^^,xo + = ^. 1.22 

n^co np[Xo) \ np[Xo) np[Xo) J 7r[u — v) 

(b) For every u,v E we have 

lim 1 L ^ w ^ ^ V \ _ ki{u) ki{v)-ki{u) ki{v) 

where Ai is the usual Airy function, and pi is the constant from ()1.18|) . 

Theorem 11.21 is similar to the main theorems Theorem 1.2 and Theorem 1.3 of It 
expresses that the local eigenvalue correlations show the universal behavior as n — oo, both 
in the bulk and at the edge, that is well-known from unitary random matrix models. So the 
result itself is not that surprising. 

To obtain Theorem 11.21 we use the Deift/Zhou steepest descent method for RH problems 
and a main tool is the three-sheeted Riemann surface associated with equation ()1.17p as 
in [Bj. There is however an important technical difference with [6]. For a > 1, the branch 
points of the Riemann surface are all real, and they correspond to the four edge points 
izi, ±2:2 of the support as described in Theorem 11.11 (c). For a < 1, two branch points 
are purely imaginary and they have no direct meaning for the problem at hand. The other 
two branch points are real and they correspond to the edge points -^zi as in Theorem 11.11 
(a). See Figure 1 for the sheet structure of the Riemann surface. The branch points on the 
non-physical sheets result in a non-trivial modification of the steepest descent method. As 
already mentioned before, one of the steps involves a global opening of lenses, and this is 
the main new technical contribution of this paper. 

The rest of the paper is devoted to the proof of the theorems with the Deift/Zhou steepest 
descent method for RH problems. It consists of a sequence of transformations which reduce 
the original RH problem to a RH problem which is normalized at infinity, and whose jump 



6 



ALEXANDER I. APTEKAREV, PAVEL M. BLEHER, AND ARNO B.J. KUIJLAARS 



matrices are uniformly close to the identity as n oo. In this paper there are four trans- 
formations Yi-*U^Th^S^R. A main role is played by the Riemann surface ()1.17p 
and certain A-functions defined on it. These are introduced in the next section, and they 
are used in Section 3 to define the first transformation Y U. This transformation has the 
effect of normalizing the RH problem at infinity, and in addition of producing "good" jump 
matrices that are amenable to subsequent analysis. However, contrary to earlier works, some 
of the jump matrices for U have entries that are exponentially growing as n — > oo. These 
exponentially growing entries disappear after the second transformation U ^ T in Section 
4 which involves the global opening of lenses. 

The remaining transformation follow the pattern of jH], JH]; Q7j and other works. The 
transformation T \—>- S in Section 5 involves a local opening of lenses which turns the remain- 
ing oscillating entries into exponentially decaying ones. Then a parametrix for S is built in 
Sections 6 and 7. In Section 6 a model RH problem is solved which provides the parametrix 
for S away from the branch points, and in Section 7 local parametrices are built around each 
of the branch points with the aid of Airy functions. Using this parametrix we define the final 
transformation S ^ R in Section 8. It leads to a RH problem for R which is of the desired 
type: normalized at infinity and jump matrices tending to the identity as n — > oo. Then R 
itself tends to the identity matrix as n — > cxd, which is then used in the final Section 9 to 
prove the Theorems 11.11 and 11.21 



2. The Riemann surface and A-functions 

We start from the cubic equation (jl.l7p which we write now with the variable z instead 
of X 

e-ze + il-a')^ + za' = 0. (2.1) 

It defines a Riemann surface that will play a central role in the proof. The inverse mapping 
is given by the rational function 

z = ^ 2.2 

There are four branch points ±zi, ±iz2 with zi > Z2 > 0, which can be found as the images 
of the critical points under the inverse mapping. The mapping ()2.2|) has three inverses, C,j{z), 
j = 1,2,3, that behave near infinity as 

U^) = z-- + Oil/z'), 

z 

U^) = a + ^ + 0{l/z'), (2.3) 
2z 

Uz) = -a+^^ + Oil/z'). 

The sheet structure of the Riemann surface is determined by the way we choose the analytical 
continuations of the ^j's. 

It may be checked that has an analytic continuation to C\ [— zi, zi], which we take as the 
first sheet. The functions .^2 and ^3 have analytic continuations to C \ \([0, zi] U [—iz2, iz2]) 
and C \ ([—2:1, 0] U [— ^^2, 1^2]), respectively, which we take to be the second and third sheets, 
respectively. So the second and third sheet are connected along [—iz2,iz2], the first sheet is 
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connected with the second sheet along [0, zi], and the first sheet is connected with the third 
sheet along [— 2;i,0], see Figure 1. 



y 



Figure 1. The Riemann surface - z^'^ + (1 - + -^^^ = 0. 



We note the jump relations 



CiT = 6± on(0,zi), 
Ci=F = 6± on (-^1,0), 
6t = 6± on {-iz2,iz2). 



(2.4) 



The A-functions are primitives of the ^-functions Xj{z) — ^j{s)ds, more precisely 

J zi 

X2{z) = f ^2{s)ds (2.5) 

J zi 

A3(^) = / Us)ds + X,_{-z,) 

The path of integration for A3 lies in C\ ((—00, 0] U [—iz2, iz2]), and it starts at the point —zi 
on the upper side of the cut. All three A-functions are defined on their respective sheets of 
the Riemann surface with an additional cut along the negative real axis. Thus Ai, A2, A3 are 
defined and analytic on C\{—oo, Zi], C\((— 00, Zi]Ll[—iz2, iz2]), and C\((— 00, 0]U[— i^2, ^^2]), 
respectively. Their behavior at infinity is 

X,{z)^^z^-\ogz + ii + 0{l/z) 



X2{z) ^ az + -\ogz + £2 + 0(1/ z) 
X^iz) = -az+^logz + £3 + 0{l/z) 



(2.6) 
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for certain constants j = 1, 2, 3. The Aj's satisfy the following jump relations 





— A2± 




on 


(0,2:1), 


Ai- 


= A3+ 




on 


(-^1,0), 


Ai+ 


= A3_ 


— vri 


on 


(-^1,0), 


A2=F 


= A3± 




on 


(0,^22), 


A2T 


= A3± 


— Tci 


on 


{-iz2,0), 


Ai+ 


= Ai_ 


- 2m 


on 


(—00, —Zi 


A2+ 


= A2- 


+ ni 


on 


(-00,0), 


A3+ 


= A3- 


+ TT'i 


on 


(—00, —Zi 



(2.7) 



where the segment {—iz2,iz2) is oriented upwards. We obtain ()2.7j) from ()2.4|1 . ()2.5|1 . and 
the values of the contour integrals around the cuts in the positive direction 



^i{s)ds = — 27ri, 
which follow from (j2.3j) . 



^2{s)ds 



^3{s)ds 



Til, 



Remark: We have chosen the segment [—^-22 5 ^2:2] as the cut that connects the branch points 
-^iz2. We made this choice because of symmetry and ease of notation, but it is not essential. 
Instead we could have taken an arbitrary smooth curve lying in the region bounded by the 
four smooth curves in Figure 2 (see the next section) that connect the points Xq, ^-225 —Xq, and 
—iz2. For any such curve, the subsequent analysis would go through without any additional 
difficulty. 



3. First transformation Y y-^U 
We define for z e C \ (M U [-iz2, ^^2]), 

U{z) = diag(e-"^\ e-"^^ e-"^«)r(z) diag(e"(^i(")"5^'), e"(^2W-a.)^ ^n{x,(z)+az)y ^g^^^^ 

This coincides with the first transformation in [Hj. Then U solves the following RH problem. 

• t/ : C \ (M U [-iz2, iz2]) C^^^ is analytic. 

• U satisfies the jumps 

(gn(Ai+-Ai_) gn{A2+-Ai_) gn(A3+-Ai_)\ 

en{A2+-A2-) Q ^^^^ (^3 2) 

gn(A3+-A3-) J 

and 

/I \ 

U+ = U^\0 e"(^2+-A2-) on[-tZ2,iZ2]. (3.3) 

\0 e"(^3+'^3-) / 

• U{z) = I + 0(1/ z) as z ^ 00. 

The asymptotic condition follows from ()1.13p . ()2.6|) and the definition of U. 



on (—00, —Zi\ 
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The jump on the real hne ()3.2j) takes on a different form on the four intervals (—00, — ^i], 
-zi, 0), (0, zi], and [zi, 00). Indeed we get from ()2.7|) . ()3.2|) . and the fact that n is even, 

gn{A2+-Ai_) gn(A3+-Ai_) 

1 

1 

^gn(Ai+-Ai_) gn{A2+-Ai_) 

1 


gra(Ai+-Ai_) 

Q gn{A2+-A2-) 



g«(A2-Ai) gn(A3-Ai) 
f/_ I 1 

,0 1 



U+ = f/_ 



f/4 



1 



gn(A3+-A3-) 
gn(A3-Ai_)^ 



1 



on [zi, 00) 



on (-2:1,0) 
on (0, zi) 



(3.4) 
(3.5) 
(3.6) 
(3.7) 



Now to see what has happened it is important to know the sign of Re (Aj — A^) for j 7^ k. 
Figure 2 shows the curves where ReAj = Re A/;. 



1 - 



E 



-2 - 





Re z 



- - ReX^ 


= ReXj 


_ , ReX^ 


= ReA.3 


ReX^ 


= ReX3 








y 




y 










y 




y 






Figure 2. Curves where ReAi = ReA2 (dashed lines), ReAi = Re A3 
(dashed-dotted lines), and ReA2 = Re A3 (solid lines). This particular fig- 
ure is for the value a = 0.4. 

From each of the branch points ±zi, ±iz2 there are three curves emanating at equal angle 
of 27r/3. We have Re Ai = Re A2 on the interval [0, zi] and on two unbounded curves from zi. 
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Similarly, Re Ai = Re A3 on the interval [—zi, 0] and on two unbounded curves from —Zi. We 
have Re A2 = Re A3 on the curves that emanate from ±iz2- That is, on the vertical half- lines 
[iz2^ +ioo) and (—200, —iz2\ and on four other curves, before they intersect the real axis. The 
points where they intersect the real axis are ±Xo for some Xq G (0,2:1). After that point we 
have Re Ai = Re A3 for the curves in the right half-plane and Re Ai = Re A2 for the curves in 
the left half-plane. Figure 2 was produced with Matlab for the value a = 0.4. The picture 
is similar for other values of a G (0, 1). As a — > 0+ or a —>■ 1 — , the imaginary branch points 
±iz2 tend to the origin. 

Using Figure 2 and the asymptotic behavior (|2.6p we can determine the ordering of Re Aj, 
j = 1,2,3 in every domain in the plane. Indeed, in the domain on the right, bounded by 
the two unbounded curves emanating from zi, we have ReAi > ReA2 > Re A3 because of 
()2.6|) . Then if we go to a neighboring domain, we pass a curve where Re Ai = Re A2, and so 
the ordering changes to Re A2 > Re Ai > Re A3. Continuing in this way, and also taking into 
account the cuts that we have for the Aj's, we find the ordering in any domain. 

Inspecting the jump matrices for U in ()3.3|) - ()3.7|) . we then find the following: 

(a) The non-zero off-diagonal entries in the jump matrices in (j3.4|) and (j3.7|) are expo- 
nentially small, and the jump matrices tend to the identity matrix as n — > 00. 

(b) The non-constant diagonal entries in the jump matrices in (|3.5j) and (|3.6p have mod- 
ulus one, and they are rapidly oscillating for large n. 

(c) The (1, 2)-entry in the jump matrix in ()3.5|) is exponentially decreasing on {—zi, — xq), 
but exponentially increasing on (— xo,0) as n — > 00. Similarly, the (l,3)-entry in 
the jump matrix in ()3.6|) is exponentially decreasing on {xo,zi), and exponentially 
increasing on (0,a;o). 

(d) The entries in the jump matrix in (j3.3|) are real. The (2, 2)-entry is exponentially 
increasing as n — > 00, and the (3, 3)-entry is exponentially decreasing. 

The exponentially increasing entries observed in items (c) and (d) are undesirable, and 
this might lead to the impression that the first transformation Y U was not the right 
thing to do. However, after the second transformation which we do in the next section, all 
exponentially increasing entries miraculously disappear. 

4. Second transformation U y-^T 

The second transformation involves the global opening of lenses already mentioned in the 
introduction. It is needed to turn the exponentially increasing entries in the jump matrices 
into exponentially decreasing ones. 

Let S be a closed curve, consisting of a part in the left half-plane from —iz2 to iz2., 
symmetric with respect to the real axis, plus its mirror image in the right half-plane. The 
part in the left half-plane lies entirely in the region where Re A2 < Re A3 and it intersects 
the negative real axis in a point —x* with x* > Zi, see Figure 3. So S avoids the region 
bounded by the curves from ±iz2 to ixo- In a neighborhood of iz2 we take S to be the 
analytic continuation of the curves where ReA2 = Re A3. As a result, this means that 

A2 — A3 is real on S in a neighborhood of iz2- (4-1) 

This will be convenient for the construction of the local parametrix in Section 7. 
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3r 





Re z 



Figure 3. Contour S which is such that ReA2 < Re A3 on the part of S in 
the left half-plane and Re A2 > Re A3 on the part of S in the right half-plane. 



The contour E encloses a bounded domain and we make the second transformation in that 
domain only. So we put T = U outside S and inside S we put 



1 
T = U [0 1 1 for Re 2 < inside E 



_e"(^2-A3) I 

(4.2) 

'10 0^ ^ ' 

T = U\0 1 -e"(^^-^2) I for Re 2 > inside E. 
1 



Then T is defined and analytic outside the contours shown in Figure 4. Using the jumps 
for U and the definition ()4.2|) we calculate the jumps for T on any part of the contour. We 
get different expressions for six real intervals, for the vertical segment [—iz2,iz2], and for 
S (oriented clockwise) in the left and right half-planes. The result is that T satisfies the 
following RH problem. 

• T : C \ (M U [-iz2, iz2] U S) ^ C^^^ is analytic. 



ALEXANDER I. APTEKAREV, PAVEL M. BLEHER, AND ARNO B.J. KUIJLAARS 




Figure 4. T has jumps on the real hne, the interval [—iz2,iz2] and on E. 



T satisfies the following jump relations on the real hne 



gn(A2+-Ai_) gn(A3+-Ai-)'« 



T+ = T_ 



T+ = T_ 



1 
^0 1 

'l e"^^^+~'^i" 
1 
,0 1 

, gn(Ai+-Ai_) g ^ 

1 

e"(^3+-^3-) 

^gn(Ai+-Ai_) 2 

gn(A2+-A2-) 

1 

'l e"(^2-Ai) o\ 
1 on[zi,x*) 

'j^ gn(A2-Ai) gn(A3-Ai) 

1 I on [zi, oo). 

.0 1 



on (— oo, — a;* 
on {—X*, —zi] 

on (-2:1,0) 
on (0, zi) 



(4. 
(4. 
(4. 
(4. 
(4. 
(4. 
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The jump on the vertical segment is 



T+ = T_ 10 1 I on[-iZ2,iz2]. (4.9) 



—1 e"(A3+-A3_) 



The jumps on E are 



T+ = T_|0 1 0| on {;z G S I Rez < 0} (4.10) 

Q gn(A2-A3) I 

'10 \ 

1 eniXs-X2)\ on {z e S I Rez > 0}. (4.11) 
.0 1 J 

• T{z) =1 + 0(1/ z) as z ^ cx). 

Now the jump matrices are nice. Because of our choice of E we have that the jump 
matrices in ()4.10|) and 1)4.111) converge to the identity matrix as n ^ oo. Also the jump 
matrices in ()4.Hj) . ()4.4j) . ()4.7)1 and ()4.8|1 converge to the identity matrix as n ^ oo. The 
(3, 3)-entry in the jump matrix in ()4.9|1 is exponentially small, so that this matrix tends to 
'l 0^ 
1 

-1 0^ 

The jump matrices in fl4.6|) and ()4.7|) have oscillatory entries on the diagonal, and they 
are turned into exponential decaying off-diagonal entries by opening a (local) lens around 
{—Zi,Zi). This is the next transformation. 

5. Third transformation T S 

We are now going to open up a lens around {zi,zi) as in Figure 5. There is no need to 
treat special point. 

The jump matrix on {—Zi, 0), see ()4.5p . has factorization 

^gn(Ai-A3)+ g 

TZ^T+ = I 1 

e"(^i-^3)- 



1 0\ / 1\ / 1 
10 10 

n(Ai-A3)- 1/1-1 0/ \e"(^i^^3)+ 1, 




(5.1) 



and the jump matrix on (0, Zi), see ()4.6|) . has factorization 



e 



n(Ai-A2) 



+ 



TZ^T+ = I 





1 0\ / 1 0\ / 1 0' 

,n{Ai-A2)- I -1 e"(^i-^2)+ I 

1/\00 1/\ 01 



(5.2) 
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Figure 5. Opening of lens around [—zi, z\\. The new matrix- valued function 
S has jumps on the real line, the interval \—iz2^iZ'^^ on S, and on the upper 
and lower lips of the lens around [— -Zi, z^. 

We open up the lens on [— zi,2i] and we make sure that it stays inside E. We assume 
that the lens is symmetric with respect to the real and imaginary axis. The point where the 
upper lip intersects the imaginary axis is called iy* . Then we define S = T outside the lens 
and 



1 0' 
S=T { 10 

_gn(Ai-A3) g ^ 
1 0^ 

S=T { 10 

e"(Ai-A3) 1^ 

1 0^ 

S = T \ -e'^(^i-^2) I 

1, 

1 0^ 

S = T \ e"(^i-^2) I 

1, 



in upper part of the lens in left half-plane, 
in lower part of the lens in left half-plane, 

in upper part of the lens in right half-plane, 
in lower part of the lens in right half-plane. 



(5.3) 



Outside the lens, the jumps for S are as those for T, while on [—zi, zi] and on the upper 
and lower lips of the lens, the jumps are according to the factorizations ()5.1|) and ()5.2|) . The 
result is that S satisfies the following RH problem 

• 5* is analytic outside the real line, the vertical segment [—iz2,iz2], the curve E, and 
the upper and lower lips of the lens around [—zi, zi]. 
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• S satisfies the following jumps on the real line 

gn(A2+-Ai_) gn(A3+-Ai_) 

S+^S^ 1 
\0 1 

(1 e'^*^'^^+~'^i-^ 
1 
1 



1' 
10 
-1 0, 



1 

5+ = 5_ ( -1 

1 

I gn(A2-Ai) g 

S_\0 1 

1 

]^ gn(A2-Ai) gn(A3-Ai) 

5_ I 1 

1 



on (— oo, —X*] 
on {—X*, —Zi] 
on (-2;i,0) 
on (0, zi) 

on [zi,x*) 



on [x*, oo). 



S has the following jumps on the segment [—iz2, iz2\, 



1 
,5+ = I 1 



—1 g"('^3+-A3-) 

1 

1 

gn(Ai-A3-) _^ gn(A3+-A3-) 

1 

1 

_gn(Ai-A3-) _^ gn(A3+-A3-) 



10 
,S+ = I 1 



—1 g'i(A3+-A3-) 



on {-iz2, -iy*) 

on {-iy*,0) 
on (0,i|/*) 
on {iy*,iz2). 



The jumps on E are 



'1 0^ 
S^^S^iO 1 

^0 gn(A2-A3) I 

'10 
5+ = 5_ I 1 e"(^3-A2) 
,0 1 



on {2; e E I Re z < 0} 
on {2; e S I Re;2 > 0}. 



(5.4) 
(5.5) 
(5.6) 
(5.7) 
(5.8) 
(5.9) 



(5.10) 
(5.11) 
(5.12) 
(5.13) 



(5.14) 
(5.15) 
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Finally, on the upper and lower lips of the lens, we find jumps 



S. 



+ 



s. 




gn(Ai-A2) I g 



gn(Ai-A3) g ]^ 



1 



1 
1 



1 




on the lips of the lens in the left half-plane 



on the lips of the lens in the right half-plane. 



(5.17) 



(5.16) 



• S{z) = I + 0{l/z) as z ^ oo. 

So now we have 14 different jump matrices ()5.4j) - ()5.17j) . As n ^ oo, all these jumps have 
limits. Most of the limits are the identity matrix, except for the jumps on {—zi, zi), see ()5.6p 
and ()5.7p . and on {—iz2,iz2), see ()5.10|) - ()5.13p . In the next section we will solve explicitly 
the limiting model RH problem. The solution to the model problem will be further used in 
the construction of parametrix away from the branch points. 



• N{z) = I + 0{l/z) as z ^ oo. 

To solve the model RH problem we lift it to the Riemann surface (12.11) with the sheet 
structure as in Figure 1, see also jB], |21| where the same technique was used. Consider to 
that end the range of the functions C,k on the complex plane, Qk = ^k{'C) for k = 1,2,3. 
Then Qi, Q2, ^3 give a partition of the complex plane into three regions, see Figure IHl In 
this figure q, p and po are such that 



6. Parametrix away from branch points 



The model RH problem is the following. Find such that 

• : C\{[-zi,zi] U [-iz2,iz2]) C^""^ is analytic. 

• satisfies the jumps 




(6.1) 



(6.2) 



(6.3) 



Q = = 6(2:1) = = -6(-Zi), 

= -6(^^2) = -^(^^2) = 6(-^^2) = 6(-i^2), 

2po = ei+(0) = -6-(0). 



(6.4) 



LARGE n LIMIT OF GAUSSIAN RANDOM MATRICES WITH EXTERNAL SOURCE, PART II 17 




-iPO 

Figure 6. Partition of the complex ^-plane. 



Let r be the boundary of Qi. Then we have 

ei±([-zi,zi]) = rn{±imz>o}, 

6-(H^2,0]) = [ip,ipo], ^2-{[0,iz2]) = [-ipo,-ip], (6.5) 
e3-(H^2,0]) = [ip,Q], e3-([0,^^2]) = [0,-ip], 

and C2±{iy) = CsTi'^v) < y < -22- According to our agreement, on the interval 

—'>'Z2 < y < "^-22 the minus side is on the right. 

We are looking for a solution N in the following form: 

fmiUz)) iVi(6(^)) N,{Uz))\ 

N{z) - iV2(6(^)) iV2(6(^)) N^iU^)) , (6.6) 
\Ns{U^)) NsiUz)) NsiUz))/ 

where Ni{^), A^2(0) -^3(0 three scalar analytic functions on C\(ru[— ipo; iPo])- To satisfy 
the jump conditions on N{z) we need the following jump relations for Nj{^), j — 1, 2, 3: 

Nj+iO = iV,_(0, e e (r n {Imz < 0}) U [-ipo, -ip] U [ip, ipo], 

Nj40 = -Nj-iO, e e (rn {imz > 0}) u [-ip,ip]. ^ ' ^ 

So the N/s are actually analytic across the curve F in the lower half-plane and on the 
segments [ip,ipo] and [—ipo,—ip]- What remains are the curve T in the upper half-plane 
and the segment [—ip,ip\, where the functions change sign. Since ^i(cxd) = 00, ^2(00) ~ a, 
^3(00) = —a, then to satisfy N{oo) — I we require 

iVi(oo) = 1, Ni{a) = 0, Ni{-a) = 0; 

Ar2(oo) = 0, N2{a) = 1, N2{-a) = 0; (6.8) 
Nsioo) = 0, Nsia) = 0, N^i-a) = 1. 
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Thus, we obtain three scalar RH problems on Ni,N2,N3. Equations ()(i.7|) - ()fi.8j) have the 
following solution: 

Ni{0= rr^S. ^2,3(0 = C2,3^==ii^==, (6.9) 



with cuts at m {Im,^ > 0} and [— ip, The constants C2,3 are determined by the equations 
A^2.3(±a) = 1. We have that 



and as in Section 6 of j6jj, we obtain C2 
problem is given by 



R{i-d) 



(6.10) 



N{Z) 



C3 = Thus, the solution to the model RH 

\ 



■_6(£)+a 



^/mMzYfi) 
„■ i'i{z)+a 



^2i?(6(^);a) v^2i?(6(2);a) 



with cuts on [— ^1,2:1] and [— ^^25 "^^2]- 



(iiz)-a 



2{z)-a 



—I 



€3(z)-a 



(6.11) 



7. Local parametrices 

Near the branch points N will not be a good approximation to S. We need a local 
analysis near each of the branch points. In a small circle around each of the branch points, 
the parametrix P should have the same jumps as S, and on the boundary of the circle P 
should match with in the sense that 

P{z) =N{z){I + 0{l/n)) (7.1) 

uniformly for z on the boundary of the circle. 

The construction of P near the real branch points ±Zi makes use of Airy functions and 
it is the same as the one given in ^ Section 7] for the case a > 1. The parametrix near 
the imaginary branch points ±iz2 is also constructed with Airy functions. We give the 
construction near iz2. We want an analytic P in a neigborhood of iz2 with jumps 

'1 0\ 
P^ = P_ \ 1 on left contour 

^10 ' 

p ( 1 e"(^3-A2) I on right contour (7.2) 




P+ = P_ 1 on vertical part. 



In addition we need the matching condition ()7.1|) . Except for the matching condition ()7.1|) . 
the problem is a 2 x 2 problem. 
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Let us consider A2 — A3 near the branch point iz2- We know that (A2 — \3){iz2) = 0, see 
fl2.7|) and since ^2 — ^3 has square root behavior at iz2 it follows that 

(A2 - A3)(^) = r - Us))ds = {z- iz^f'^z) 

J iz2 

with an analytic function h with h{iz2) 7^ 0. So we can take a 2/3-power and obtain a 
conformal map. To be precise, we note that 



arg((A2 - A3)(2?/)) = 7r/2, for y > 2:2, 



and so we define 



such that 



-(A2-A3)(^) 



2/3 



aTgf{z) = n/3, hTz = iy,y>Z2. 



(7.3) 



Then s = f{z) is a conformal map, which maps [0,'i2;2] into the ray args = — and which 
maps the parts of S near iz2 in the right and left half-planes into the rays args = and 
args = respectively. [Recall that A2 — A3 is real on these contours, see (|4.ip .] 
We choose P of the form 



'1 

P{z) = E{z)^ {n'^/^f{z)) ( e3"(^2-A3) 



.0 







e 2 








(7.4) 



where E is analytic. In order to satisfy the jump conditions ()7.2|) we want that $ is defined 
and analytic in the 
there it has jumps 



and analytic in the complex s-plane cut along the three rays args = k^, k = —1, 0, 1, and 




for arg s = 27r/3, 

for args = — 27r/3, 
for args = 0. 



(7.5) 
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Put yo{s) = Ai(s), yi{s) = uAi^us), y2{s) = u'^Ai^u'^s) with u = 27r/3 and Ai the 
standard Airy function. Then we take $ as 

/I \ 

= yo -1/2 for < args < 27r/3, 

\o y'o -2/2/ 
/I o\ 

$ = ?/o 1/1 for - 27r/3 < args < 0, (7.6) 




20F 



$= -yi -y2 \ for 27r/3 < args < 47r/3. 

yo -y'l -y'2/ 

This $ satisfies the jumps fl7.5|l . In order to achieve the matching ()7.H1 we define the prefactor 
E as 

E = NL-^ (7.7) 

with 

'1 \ /l 

1 i I (7.8) 

^0 n^/^f^/y \0 -1 

where /^/^ has a branch cut along the vertical segment [0, iz2\ and it is real and positive 
where / is real and positive. The matching condition ()7.1|) now follows from the asymptotics 
of the Airy function and its derivative 

Ai(s) = ^s-VVi^^^^ (l + 0(s-^/^)), 
Ai'(s) = - J^s^Vi^^^^ (l + 0(s-^/^)), 
as s ^ 00, I args I < vr. On the cut we have f^/^ = if]/^. Then fl7.8|) gives 




which is the same jump as satisfied by A^, see (j6.3p . This implies that E = NL^^ is analytic 
in a punctured neighborhood of iz2- Since the entries of A^ and L have at most fourth- 
root singularities, the isolated singularity is removable, and E is analytic. It follows that P 
defined by ()7.4|) does indeed satisfy the jumps ()7.2|) and the matching condition ()7.1|) . 
A similar construction gives the parametrix in the neighborhood of —iz2- 

8. Fourth transformation S ^ R 
Having constructed A^ and P, we define the final transformation by 

R{z) = S{z)N{z)~^ away from the branch points, 
R{z) = S{z)P(z)~^ near the branch points. 
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Since jumps of S and coincide on the interval {—zi, zi) and the jumps of S and P coincide 
inside the disks around the branch points, we obtain that R is analytic outside a system of 
contours as shown in Figure 7. 




Figure 7. R has jumps on this system of contours. 

On the circles around the branch points there is a jump 

R+ = R^I + 0{l/n)), (8.2) 

which follows from the matching condition ()7.H) . On the remaining contours, the jump is 

i?+ = i?_(/ + 0(e~^")) (8.3) 

for some c > 0. Since we also have the asymptotic condition R{z) = I + 0(1/ z) as z ^ oo, 
we may conclude as in Section 8] that 



R{z) 



I + O 



n{\z\ + 1) 

uniformly for z G C, see also [IHl CEl IT7| 1^. 



as n — > oo 



9. Proof of Theorems 11.11 and 11.21 

We follow the expression for the kernel as we make the transformations Y U t-^ T i-^ 
S. From p.l4|) and the transformation 1)3.11) it follows that Kn has the following expression 
in terms of U, for any a;, y G M, 



gin(a;2-y2) 

2ni{x — y) 



g-nAl+(a;)^ 





(9.1) 
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Then from ()4.2j) we obtain for y > inside the contour S, and for any x G 

^n(a:,y) = |jT^^-^(0 e-'^+^y^ O) T~\y)T+{x) i |, (9.2) 

and from ()5.3|) . we have when x,y E [0, Zi), 

KJx,y) = (-e"^i+(^) e"^2+(J') O) g-^^^+W | . (9.3) 

27r2(a; - ' I 

Since Ai+ and A2+ are each others complex conjugates on [0,2;i), we can rewrite (j9.3|) for 
x,y E [0, zi) as 



^n{h{y)-h{x)) 



e 



KJx,y) = ^ (_e"^i"i^i+fe) e-^^'^^i+^s') O) e"i°>^i+(^') , (9.4) 

2m[x — y) ' ' 

where 

h{x) = Re Xi+{x) - ^x^ . (9.5) 

Note that (j9.4j) is exactly the same as equation (5.14) in 6j. Therefore we can almost literally 
follow the proofs in Section 9 of fB| to complete the proof of Theorem 11.11 and 11.21 

Indeed as in jH] the limiting mean density ()1.15|) follows from ()9.4|) and ()8.4|) in case x > 0, 
where 

p{x) = -Im^i+(x), xeR. (9.6) 

The case x < follows in the same way and also by symmetry. Recalling that the choice of 
the cut [—iz2,iz2] was arbitrary as remarked at the end of section 2, we note that we might 
as well have done the asymptotic analysis on a contour that does not pass through 0, so that 
we obtain p.l5j) for x = as well. The statement in part (a) on the behavior of p follows 
immediately from ()9.6p and the properties of ^1 as the inverse mapping of ()2.2p . 
This completes the proof of Theorem 11.11 

The proof of part (a) of Theorem 11.21 for the case Xq > follows from (j9.4p and (j8.4j) 
exactly as in Section 9 of jB]. The case Xq < follows by symmetry, and the case Xq = 
follows as well, since we might have done the asymptotic analysis on a cut different from 
[—iz2, iz2\i as just noted above. The proof of part (b) follows as in as well. Note however 
that the proof of part (b) relies on the local parametrix at the branch point 21, which we 
have not specified explicitly in Section 7. However, the formulas are the same as those in jH] 
and the proof can be copied. This completes the proof of Theorem 11.21 
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